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Abstract 

A suitable derivative of Einstein's equations in the framework of 
the teleparallel equivalent of general relativity (TEGR) yields a conti- 
nuity equation for the gravitational energy-momentum. In particular, 
the time derivative of the total gravitational energy is given by the 
sum of the total fluxes of gravitational and matter fields energy. We 
carry out a detailed analysis of the continuity equation in the context 
of Bondi and Vaidya's metrics. In the former space-time the flux of 
gravitational energy is given by the well known expression in terms of 
the square of the news function. It is known that the energy definition 
in the realm of the TEGR yields the ADM (Arnowitt-Deser-Misner) 
energy for appropriate boundary conditions. Here we show that the 
same energy definition also describes the Bondi energy. The analy- 
sis of the continuity equation in Vaidya's space-time shows that the 
variation of the total gravitational energy is determined by the energy 
flux of matter only. 
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1 Introduction 



In spite of its wide acceptance as a theory for the dynamics of the space- 
time and of the gravitational field, Einstein's theory of general relativity, in 
the metrical formulation, does not allow the emergence of an unequivocal 
definition of gravitational energy. The definition of such quantity is surely 
important for a comprehensive understanding of the theory. In the past 
the major attempts at defining the gravitational energy were carried out by 
means of pseudotensors [1, 2, 3]. Nowadays the concept of quasilocal energy 
(i.e., the energy associated to a closed spacelike 2-surface) is being thoroughly 
investigated in the context of Hamiltonian formulations (see, for instance, 
Ref. [4] and references therein). Although this approach is promising, there 
is no ultimate consensus about it. 

Two well established notions of gravitational energy are the total energy 
of the space-time, known as the ADM (Arnowitt-Deser-Misner) energy [5], 
evaluated on a spacelike slice of a space-time that is asymptotically flat, 
and the Bondi energy [6] that describes the mass of radiating systems in 
asymptotically flat space-times. The latter is usually defined at null infinity. 
A remarkable feature of the Bondi energy is that it is related (by means of 
the field equations) to the loss of the total mass of a source that radiates 
gravitational waves. The loss of total mass is given in terms of minus the 
square of the news function, and therefore the mass of the source can only 
decrease. 

A definition for the gravitational energy-momentum has been investigated 
in the realm of the teleparallel equivalent of general relativity (TEGR) [7, 8] 
(Ref. [9] considers the TEGR as a particular case of metric-affine theories of 
gravity, and provides a large number of references on the subject). Attempts 
to define the gravitational energy in the context of teleparallel theories of 
gravity were first considered by M0ller [10], who noticed that an expression 
for the gravitational energy density could be constructed out of the torsion 
tensor. Arguments in favour of the localizability of the gravitational energy 
are presented in Refs. [7, 8, 11], where it is argued that it is possible to 
identify the gravitational energy density with an appropriate scalar density, 
defined in terms of a set of global tetrad fields, which is a total divergence. 
Therefore the gravitational energy density in consideration transforms as a 
scalar density under coordinate transformations. We recall that the local- 
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izability of the gravitational energy was already supported by Bondi [12], 
who analyzed the possible forms of energy transfer of the gravitational field. 
He claimed that in relativity a nonlocalizable form of energy is inadmissible, 
because all forms of energy contribute to gravity, and therefore its location 
can in principle be found. 

However, it is reasonable to expect the localization of the gravitational 
energy to depend on the reference frame [11]. The Principle of Equivalence 
is based on the equality of inertial and gravitational masses. Alternatively, 
an accelerated frame can be locally considered as a rest frame with the addi- 
tion of a certain gravitational field. Therefere our perception of the strength 
of the gravitational field on nearby bodies depends clearly on our reference 
frame, and so does the gravitational energy as measured on the same frame. 
Reference frames are described in terms of tetrad fields [13, 14], which are the 
basic field quantities of the TEGR. The total energy of a space-time, how- 
ever, is not expected to depend on a particular frame (the total gravitational 
energy of any asymptotically flat space-time, in the TEGR, is precisely the 
ADM energy [7]). The dependence of the localization of the gravitational 
energy on the frame is a manifestation of the fact that gravitational and 
inertial forces are of the same nature. 

In this article we address the issue of the variation of the gravitational 
energy and the corresponding total flux of energy at spacelike infinity. The 
success of the present analysis may qualify the TEGR as an adequate frame- 
work to investigate more intricate astrophysical models of energy flux and 
gravitational waves emission. We conclude that the energy-momentum def- 
inition in consideration here provides a unique framework to analyze both 
the ADM and the Bondi energy. We will show that for the Bondi space-time 
metric tensor the decrease in time of the total gravitational energy causes the 
flux of gravitational energy only (by means of the emission of gravitational 
waves), whereas in Vaidya's space-time [15, 16] the loss of gravitational en- 
ergy is related to a pure flux of matter energy. The variation in time of the 
gravitational energy is determined by the sum of the gravitational and matter 
energy fluxes. From the point of view of Bondi [12], these fluxes characterize 
the intangible and tangible forms of gravitational energy transfer, respec- 
tively (the latter being determined by the energy- momentum tensor). 

In section II of the article we recall the simple derivation of the continuity 
equation for the gravitational energy-momentum. In section III we present 
the construction of the simplest set of tetrad fields, defined in terms of the 
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retarded time u — t — r, that describes the Bondi space-time. The Bondi 
energy, the loss of total mass and the total flux of gravitational energy at 
spacelike infinity are obtained in section IV. In section V we consider the 
Vaidya space-time metric tensor and relate the loss of the total mass of the 
source to a pure matter energy flux. Finally in section VI we present our 
conclusions. 

Notation: space-time indices //, u, ... and S0(3,l) indices a, b, ... run from 
to 3. Time and space indices are indicated according to /i — 0, i, a = (0), (i). 
The tetrad field e a M yields the definition of the torsion tensor: T a ^ = 
dfj,e a v — d u e a M . The flat, Minkowski space-time metric is fixed by r\ ab = 
e a ^e bu g^ = (- + ++). 



2 The fluxes of gravitational and matter fields 
energy 

The Lagrangian density for the graviational field in the TEGR, in the pres- 
ence of matter fields, is given by 



L(e a ,) = -ke{-T abc T abc + l -T abc T bac -T a T a )-L m 

= -keZ abc T abc -L m , (1) 

where k = l/(167r) and e = det(e a At ). The tensor S a6c is defined by 

Y^abc ^^rpabc _|_ rpbac rpcab^ _|_ _^jy ac J 1 ' > jf^T '} (2) 

4 2 
and T a = T b b a . The quadratic combination T, abc T abc is proportional to the 
scalar curvature R(e), except for a total divergence [17]. L m represents the 
Lagrangian density for matter fields. The field equations for the tetrad field 
read 

e aX e b MeZ bX ») - e(S^ a T bufl - -^T,^) = j^eT afl , (3) 

where 
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It is possible to prove by explicit calculations that the left hand side of Eq. 
(3) is exactly given by \ e \R a[l {e) - \e a[l R{e)\. 

Multiplication of Eq. (3) by the appropriate inverse tetrad fields yields 
the latter equation in the form 



d u {-AkeL aX ») = -kee^i^T^ - 5fi bdc T bcd ) - ee a ^ . (4) 

As usual, tetrad fields convert space-time into Lorentz indices and vice- versa. 
By making A = j (i.e., by restricting the space-time index A to assume only 
space values) we obtain 



-do(-4keY?°*) - d k {-Akell ak i) = - kee^^T^ - 5^ bcd T bcd ) 

- ee\T^. (5) 

Now we take the derivative of the equation above with respect to j. The 
second term on the left hand side vanishes because of the antisymmetry 
property S 61 ^ = — S ajA: . The resulting equation reads 



-dodji-AkeE^) = - ^[ee^(4S^T 6cM -^S 6cd T 6cd )] 

- dj(ee\T^). (6) 

We recall that in the Hamiltonian formulation of the TEGR [18] the 
momentum canonically conjugated to the tetrad components e a j is given by 
n a ^ = —AkeTj a0 ^ and that the gravitational energy-momentum P a contained 
within a volume V of the three-dimensional spacelike hypersurface is defined 

by [7] 

P a = -J ^xdjU aj . (7) 

If no condition is imposed on the tetrad fields, P a transforms as a vector 
under the global SO(3,l) group. (The Hamiltonian formulation of the TEGR 
has also been extensively discussed by Blagojevic and Nikolic [19].) 



4 



By integrating Eq. (6) on a volume V of the three-dimensional space we 
arrive at 



d_ 

dt 

where 



d 3 xd i U aj 

v 



- , (8) 



* a g = k[ dS 3 [ee^(4^T bctl - 5fc bcd T bcd )} , (9) 
J s 

is a component of the gravitational energy-momentum flux, and 



a m = f g dS j (ef ll T*>), (10) 

is the a component of the matter energy-momentum flux. S represents the 
spatial boundary of the volume V. Therefore the loss of gravitational energy 
p(°) = E is governed by the equation 

*E = _$(o) _ $(0) (n) 
dt 9 m ' K ' 

Equation (8) has been derived previously in the Lagrangian formulation in 
Ref. [11], without the addition of matter fields. It does not depend on details 
of any Hamiltonian formulation. The time variable t may also stand for the 
retarded time u. We note that in the analysis of the Hamiltonian formulation 
of the gravitational field in null surfaces, in the framework of the TEGR [20], 
the energy definition given by Eq. (7) was already anticipated. However 
the investigation of the gravitational energy for the Bondi space-time metric 
tensor was first carried out by imposing the time gauge condition [21] (use 
was not made of the retarded time), a framework that is not quite suitable 
for such analysis. In the following sections we will show that $^ plays a 
major role to Bondi's radiating space-time, since it yields the well known 
expression for the gravitational energy loss, whereas only $^ is relevant to 
Vaidya's space-time. 



3 The tetrad field for the Bondi space-time 

The tetrad field for the flat space-time, in cartesian coordinates, is given by 
e a n(ct,x,y, z) = 5". By making a coordinate transformation t — > u, where 
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u = t — r/c is the retarded time and r = \Jx 2 + y 2 + z 2 , we obtain 



e a v(u,x,y,z) 



[1 



Vo 



x/rc 
1 





y/rc 

1 




z/rc\ 



1 J 



(12) 



where a and fi label the line and column indices, respectively. In spherical 
coordinates we have 



1/c \ 

sin 9 cos r cos cos — rsin#sin0 , , 

sin 9 sin r cos 6? sin r sin # cos 

cos — r sin 9 / 

In Eqs. (12) and (13) only c stands for the velocity of light. 

The set of tetrad fields given by Eq. (12) defines a holonomic transforma- 
tion dq a = e a ^dx^ between the coordinates — u + r/c, = x, q^ = y 
and = z that describe a reference space-time [7], and the space-time coor- 
dinates (u, x, y, z). In the form given by Eq. (12) the tetrad fields violate the 
time gauge condition because ^ 0, and therefore eu\ 7^ 0. A nontrivial 
manifestation of the gravitational field takes place when the transformation 
dq a = e a ^dx^- becomes anholonomic [7] . 

It is known that the use of the retarded time u clearly simplifies the de- 
scription of the Bondi space-time, since in this case the null surface condition 
g 00 = can be easily achieved. However it is not possible to impose at the 
same time the condition g 00 = and the time gauge condition. 

In the flat space-time the violation of the time gauge condition may be 
associated to a boost of the physical space-time with coordinates x^ with 
respect to the reference space-time with coordinates q a [7]. The boost singles 
out a particular direction in space. If the time gauge condition is to be 
violated in the description of a radiating gravitational field, then it should not 
privilege any direction in space, and ideally the violation should be isotropic. 
This feature is manifest in the flat space-time tetrad field given by Eq. (12), 
and it will characterize the simplest set of tetrad fields that describes the 
Bondi space-time metric tensor, as we will see. 

The Bondi metric tensor is given by 



1 M (u, r, 9, 0) 



/l 




Vo 
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ds 2 = - {Y- e ^-U 2 r 2 e 2 ^du 2 -2e^dudr 

- 2Ur 2 e 2 ^ dud9 + r 2 [e^d9 2 + e~ 27 sin 2 9 d^, (14) 

where (r, 9, 0) are the usual spherical coordinates and u is the retarded time. 
The functions in the equation above have the following asymptotic behaviour 



4r 2 



7 
V 
r 



c 

- + 

r 



2M 



dd 



+ d cos 9 



fdc\- 
\d9j 

+ ■■■ 



4e(|)oat« 



U 



--c 2 (1 + 8 cot 2 9 
1 fdc 

-Am +2ccote 

1 / dc 
+ —[2d + 3c— cot 9 + 4c 2 cot 9 ) 



(15) 



where M = M(u, 9) and d = d(u, 9) are the mass aspect and the dipole 
aspect, respectively, and from the function c(u, 9) we define the news function 
c = (dc(u, 9))/du. The inverse metric tensor is given by 



/ 



«f>,r,0,0) 



-2/3 



p -2(3V 





p- 2 i 







e 2 ^ 



(16) 



-e-^U e - 

i — 

V ^f-r-J 

\ r A sin ' 

and the determinant of e a M reads e = y/—g = e 2/3 r 2 sin 2 9. 

The simplest set of tetrad fields e a ^(u, r, 9, <fi) that leads to Eq. (14) 
is precisely the one for which the violation of the time gauge condition is 
isotropic. It reads 



7 



e, 



/ -A -B 

C cos 9 cos B sin cos re 7 cos cos 

C cos # sin B sin sin re 7 cos # sin 

\ — Csin# Bcos9 —re 1 sin 9 



-re 7 sin 9 sin 
re~ 7 sin# cos0 











(17) 



The functions A, B and C are defined by 



B 



C 



A 




2 



(18) 



It is possible to show by means of a coordinate transformation that 



Therefore the feature above mentioned regarding the violation of the time 
gauge condition is manifest in Eq. (19) (of course the functions (3 and V 
depend on 9). In this respect Eq. (19) shares a similarity with Eq. (12). 
Moreover the spatial components of e afl in Eq. (17), in cartesian coordinates, 
satisfy 



For the flat space-time tetrad field Eq. (20) ensures that the physical space- 
time with coordinates is not rotating with respect to the reference space- 
time with coordinates q a . This condition is required to hold for a set of tetrad 
fields that describe an arbitrary gravitational field [7]. 




(19) 



e {i)j (u,x,y,z) =e {j)i {u,x,y,z) . 



(20) 



8 



The set of tetrad fields given by Eq. (17) is fixed precisely by the six 
conditions determined by Eqs. (19) and (20). We propose that the isotropic 
violation of the time gauge condition, given in the present context by Eq. 
(19), together with the condition defined by Eq. (20) fix the structure of 
any set of tetrad fields for radiating space-time metric tensors determined 
in terms of the retarded time u, in which case the imposition of the time 
gauge condition is troublesome. We argue that such tetrad field structure is 
adapted to static observers at spacelike infinity (because of Eq. (20)) and 
is suitable to investigate the properties of the energy flux of a system that 
emits gravitational waves, as we will see. 

Before closing this section we note that the idea of a preferred frame has 
been recently considered in connection with investigations either in cosmol- 
ogy or in attempts to the quantum theory of gravity [22]. The motivation 
behind such attempts to a preferred frame are quite different from our pur- 
poses, and even the approach itself is conceptually and mathematically differ- 
ent from the what has been presented above. Nevertheless it is worth noting 
that similar ideas emerge in several distinct analyses in general relativity. 



4 Bondi energy and the energy loss 

By considering a = (0) in Eq. 7 we have an expression for the energy of 
the gravitational field contained within a volume V of the three-dimensional 
spacelike hypersurface, bounded by a surface S of constant radius r, 



The tensor S^ 01 = S 001 + e^jff 01 is evaluated by means of Eq. (2). 
It simplifies to 



E = — 



[ d 3 xd,jU^ 



JV^oo 




(21) 
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E (o)oi = I e (o) ( T ooi + ^i TO) 

+L ( -°\(T 101 + g llr t° - g^T 1 ) . (22) 

We need to calculate all components T\^ u = e a \T afll/ of the torsion tensor 
obtained from Eq. (17). For later purposes, we also provide the asymptotic 
behavious of these componentes with respect to the radial variable r, in the 
limit r — > oo. The nonvanishing components are given by 

Tooi = Wa'-C^-AOoB ~ O(-) 

T 10 i = BdxA ~ O(l) 

T 201 = -re^C ~ O(-) 

r 

T 002 = ^ 2 (A 2 -C 2 ) + C(9 o7 )re^ O(-) 

T 102 = Bd 2 A + BC ~ O(-) 

r 

T 202 = re 7 [(<9 o7 )re 7 -<9 2 C] ~ O(r) 

T 303 = -(<9 o7 )rV 27 sin 2 # - re~ 7 C sin cos ~ O(r) 

T 012 = ASb J B + C[e' r + (9 1 7)re 7 - J B] ~ O(-) 

T 212 = re 7 [e 7 + (d l7 )re 7 - B] ~ O(r) 

T 313 = re" 7 ^ 7 - (9i 7 )re- 7 - £] ~ 0(r) 

T 323 = -r 2 (<9 27 )e~ 27 sin 2 #-r 2 (l-e~ 27 ) sin 6 cos 9 ~ 0(r 2 ) . 

(23) 

We also need the traces T^. They read 



T° = sV^ioi 
T 1 = 



A/^212 + <? 33 T 313 ) 



9 J 001 



gWon-T. 



—9 9 J- 202 — 9 9 1 303 — 9 9 J- 323 + 9 9 L '< 



,01„33 



201; 
12 33^ 



9 9 1 212 
,l2A2r 



g'V'T, 



313 



212 
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T 2 — „01_12 T 12 22 T 12 33 T 

J — <? 5 J 101— 9 9 +212— g 9 J- 313 

T 3 = 

After the substitution of Eqs. (23) and (24) into Eq. (22) we obtain 



(24) 



£ {0)01 = -\e {0 \[g m g m g 22 T 212 + g m g ^ g ^ 



313 



+ 2 e n9 9 9 1 202 + 9 9 9 + 303 + 9 9 9 J323J • {^) 

With the help of the equation above we arrive at an expression for the total 
gravitational energy of the space-time, 



E 



— lim / 



,1 



+r 2 e~ 13 



d6—< r sin 6 
2 1 

V jd (U sin 6) 



e < + e~ 7 - 2e 



V r 



7 (7 



r / 09 

= - lim / doirsinO e 7 + e 
4 r ^°° Jo I 

Considering the last integral in Eq. (26), we note that 



(26) 



lim / d9r e 



- 2 e-P(-) *-^(Usin9) = lim d6r 2 ^-(U sin 6) . (27) 
r^oo jo \r J 06 r ^°° Jo a# 

Since the function U{9) sin6 is required to vanish at 9 = 0, n [6], we conclude 
that the integral above vanishes. Moreover we have 



e 7 + e~ 7 - 2e~ p - 



V\h 2M c 2 



in the limit r — > 00. Therefore in such limit we obtain 



E = \ d6 sind M(u,9) , 
2 Jo 



(28) 



(29) 
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which is precisely the expression of the Bondi energy. 

It is possible to show by explicit calculations that the Bondi energy is 
invariant under local Lorentz transformations, e a = A a b (x)e b ^(x), whose 
asymptotic behaviour is given by 



A (0) (o) = 1+ V 0) ( o)(l/r) 
A(°) (i) = V°) (i) (l/r) 

AW -) = *g> + D W W W + V^l/r), (30) 

where and 1 cj a b are infinitesimal quantities such that = 

— °c<j(j)(j), 1 uj a b = — 1 Ub a , and are constants. This result is no sur- 

prise since Eq. (7) can be evaluated as a surface integral, and all transfor- 
mations that preserve the asymptotic structure of the field quantities are 
allowed. Nevertheless it is possible to show by explicit calculations that in 
the limit r — > oo we have 



^212 


— T 2 12 


T313 


= 2313 


^202 


— T202 


^303 


= ^303 


T323 


= ^323 



(31) 
and therefore 



S (o)01 = -i(A (0) (0)e W + A<°> (i)e W )k° W 2 T 2 i 2 + ^W^d 

+i(A® ( o)e (0) 1 + A<°> W e« OfoVV 2 ^ + sVV^aos 
VVV^m , (32) 

from what we conclude that Eq. (21) is invariant under the transforma- 
tions given by Eq. (30). We note that in the expression above the tetrad 
components 1 behave as r° when r — > 00. 
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The flux of gravitational energy is determined by the a = (0) component 
of Eq. (9), and by making j — 1, i.e., by integrating over a surface S of 
constant radius r, and requiring r — > oo. It reads 

$f = fc / ^[ee^(4E^T bc/1 - <%E bcd T bcd )] 




+<7 U e ( > ^ ! + 5 12 e(°) ^ 2 ] , (33) 

where 



fi 1 


i 


= 4E Act1 T A(j1 - 







= 4E Act1 T Act0 , 


n 1 


2 


= 4E Act1 T A(t2 . 



The expression of the quantities above are given in the Appendix. In 
order to calculate the total flux in the limit r^oowe take into account the 
asymptotic behaviour of the following field quantities, 



9 01 - -1, 

9 U - 1, 

g 12 - 0(l/r 2 ), 

e (0) o - I, 

e (0) o - 1, 

e — > r 2 sin# . (35) 

It is easy to verify that in the limit r — > oo the first and third term on the 
right hand side of Eq. (33) cancel each other, and that the last term behaves 
as g l2 e^ 2 = 0(l/r 4 ). Therefore $(°) reduces to 

$ f = T^/ d6d(f>e{-n\) . (36) 

a 107T JS^oo 

In order to simplify the evaluation of the integral above we only consider 
terms in Eq. (A. 2) that fall off as 1/r 2 . For this purpose we take into account 
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Eqs. (23) and (35). We conclude that five terms may possibly contribute to 
the surface integral when r — > oo. By integrating in we find that in this 
limit Eq. (36) simplifies to 



By explicitly calculating the asymptotic behaviour of all quantities under in- 
tegration, and after carrying out several cancellations, we conclude that only 
one term arising from the product T202T303 in Eq. (37) yields a nonvanishing 
value to <J>(°). We obtain 



which is the well known value of the loss of mass in Bondi's space-time [6]. 
We note that no multiplicative factors had to be adjusted in Eqs. (9) or (33), 
in order to arrive at the expression above. To our knowledge it is the first 
time that the loss of mass has been obtained by means of a flux equation 
for the gravitational energy. Finally we note that since Eq.(14) is a vacuum 
solution of Einstein's equations, then $^ = and thus Eq. (11) is naturally 
verified. 



5 The Vaidya metric tensor 

The Vaidya metric [15] describes the exterior geometry of a spherically sym- 
metric star when the energy flux of matter is included. The star is consid- 
ered to be throwing away mass. A spherically symmetric matter distribution 
does not emit gravitational radiation by means of gravitational waves, as in 
Bondi's space-time. Therefore in the present context the energy flux is a 
purely matter flux. The metric is determined by the mass function m(u), 
which is a nonincreasing function of the retarded time. It reads 




(37) 




(38) 



ds 2 = -(1 - 



2m{u) 



)du 2 - 2du dr + r 2 (d6 2 + sin 9 2 dcf 2 ) . 



(39) 



r 
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If we require the function m to be constant, then the metric tensor above 
reduces to the Schwarzschild solution. 

The set of tetrad fields that yields Eq. (39) is obtained by following the 
same steps that led to Eq. (17). We require e a M to satisfy a condition similar 
to Eq. (19), i.e., the isotropic violation of the time gauge condition, together 
with Eq. (20). The tetrad field for which both conditions are satisfied is 
given by 

I -a -a' 1 \ 

a" 1 sin 6 cos <fi rcos6>cos0 — rsin6>sin0 , , 
a -1 sin sin rcos#sin0 rsin#cos0 ' 
\ a -1 cos# — rsin6> / 

2m/r) 1 / 2 . The nonvanishing components of T^ v \ = e a ^T au \ 



e aii (u,r,6,(j>) = 

where a — (1 — 
are 



m m 

1 ooi — 5 — I — ~ 

olt r z 
m 



a 2 r 2 



T 2 i2 = r(l - a v ) , 

T 313 = r(l -aT 1 ) sin 2 , (41) 

with m = dm/du. The calculations in this context are much simpler than in 
the preceeding case. We have 



E (o)oi = e (o) oS ooi + e (o) iS 



101 



- e (°) (T 001 + 9 01 T°) 

+ lg(0) i(T l0l + glljO _ g 0l T l^ (42) 



After a number of simplifications we obtain 



S (0)01 = -ie(°) (^ 2 T 212 + ^T 31 3) 



= l -{l-a). (43) 
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Therefore by considering a finite volume V, enclosed by a spherical surface 
S of radius r, we have 



E(u) = - J (PxdjllW 



= [ d9d<P(4ke^ 0)01 ) 
J s 



rll-Jl-^1. (44) 



The total gravitational energy of the space-time is obtained by making r — > 
oo. It is given by E to tai = m(u). From Eq. (44) it follows that 

dE(u) = m(u) 

du ~ ~ ^ _ 2mW ' 

Now we proceed to evaluate the fluxes $^ and The former is given 

by 

$f = k I dS 3 \ee^{AT!^T bc ,-5^ bcd T bcd )} 
J s 

= kf dS 1 e[g 01 eM Q Qi 1 + ff oi e (o) q + g ii e (o) ±Q i ±] (46) 
J s 

where S is a surface of constant radius r and the quantities are defined by 
Eq. (34). For the Vaidya metric tensor we have g 00 = 0, g 01 = — 1, g 11 = a 2 
and g 12 = 0. By taking into account the nonvanishing components of T^ uX it 
is possible to conclude that in the present case we find 

5^0 = 0, 

n\ = (47) 

Moreover we have g 01 e^ o = — ot and g ll e^ i — a, from what we conclude 
that <J>(°) = 0, as expected. 

In order to calculate the energy flux of matter we need the energy- 
momentum tensor associated to Eq. (39). It reads [16] 
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For a volume V enclosed by a surface of constant radius S we find 

a® = / ^(ee^T^) 
= / dS^ ee (0) 

= f d9 d(f){r 2 sin 0) e (°> iT 11 . (49) 
Given that T 11 = — m/(47rr 2 ), we obtain 



*<? = -- , (50) 



A 3 



2m(u) ' 



Therefore Eqs. (45) and (50) are in agreement with Eq. (11). We conclude 
that in the Vaidya space-time the energy flux of matter accounts for the 
loss of gravitational energy. By taking the limit r — > oo in Eqs. (45) and 
(50) we recover the well known result of the total flux (total luminosity L^) 
measured by an observer at infinity, = —dm/du [16]. 



6 Concluding remarks 

The flux of the energy-momentum of gravitational and matter fields has not 
been investigated so far in the literature in a unique geometrical setting. 
In particular, to our knowledge the flux of matter fields $^ defined by Eq. 
(10) has not been previously considered. The total flux of matter fields does 
not depend on the particular structure of the tetrad fields, but only on the 
asymptotic behaviour of the latter. The remarkable feature of the continuity 
equation (8) is that it is a natural consequence of Einstein's equations written 
in the TEGR. The detailed analysis of Eq. (11) in the context of Bondi and 
Vaidya space-times demonstrates the consistency of the present framework 
for investigating the energy-momentum properties of the gravitational field. 
One may then consider the application of the techniques presented here to 



17 



concrete astrophysical problems, in which case one may assume the surface 
S in Eq. (9) to be an open surface, as considered in Ref. [11]. In the 
latter reference we evaluated the energy flux of plane gravitational waves. 
It is likely that for a given set of tetrad fields a nonvanishing value of <&(°) 
indicates the emission of gravitational waves. However, Bondi argues that 
the gravitational energy can be transferred by means of two distinct ways, 
radiative transfer and inductive transfer [12]. Therefore it is possible that 
a nonvanishing value of $^ also indicates a form of inductive transfer of 
gravitational energy. This issue will be further analyzed. 

Restricting considerations to spherical surfaces S of constant radius, it is 
possible to show by explicit calculations that for the Bondi space-time the 
three momenta PW given by Eq. (7) vanish. However, we note that P^ 3 ' 
vanishes only by placing the surface S at spacelike infinity. For a spherical 
surface S of finite radius r, is nonvanishing. The physical meaning of a 
nonvanishing pW is under investigation. 

It is an interesting achievement of the present approach to the energy- 
momentum of the gravitational field that the definition given by Eq. (7) 
allows a unfied treatment of the Bondi and ADM energies. For this purpose 
we have to fix suitable expressions for the tetrad field in each case. Although 
the localization of the gravitational energy depends on the choice of the tetrad 
field, total quantities are in fact invariant under a wide class of local Lorentz 
transformations. 
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APPENDIX 

For the Bondi space-time metric tensor the quantities Q 1 * v are given by 

i — g 01 g 01 g 22 (^-Toi2Toi2 + T012P201 - -^TwiT^i 



~ z g g g - f ooi J 3i3 + 47 g g +012+: 



212 
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+2^ 1 ( 7 1 V 3 (Toi2T 31 3 - T 201 T 313 ) 
+g 12 g 12 g 22 T 212 T 212 + 2g 12 g 12 g 33 T 212 T 313 
■ig v g n g"W: U , - / W^T^ , (Al) 

o — g 01 g 01 g 22 ^3T o2^io2 + ^002^012 + ^002^201 

-2T 202 T 00 i) - 2g m g 01 g 33 T 001 T 303 + 2g ll g l2 g 22 T W2 T 212 

+2g m g 12 g 22 ^002^212 + t 2 q 2 Tiq 2 + t 2 o 2 Toi 2 ^J 
+2g m g 12 g 33 ^002^313 - Ti 02 t 303 + T012T303 
~~ ^303^201^ + ^g m g 11 g 22 {tiq 2 Tqi 2 + t 10 2T 10 2^ 

—Ti 02 T 323 + T 202 T 3 i 3 + T303T212 
-2g n g 22 g 33 (t 202 T 313 + T 21 2T3o 3 ) - 4<? 01 g 22 l g 33 T 202 T 303 , ( A.2) 

2 = g 01 g 01 g 22 (^—2Too 2 T 2 i 2 — t 2 qiT 2 q 2 + t 2 q 2 T\o 2 
+^202^012^ + 2g m g 01 g 33 ^102^303 — ^001^323 - 2002^313^ 
+2g 01 g 12 g 33 t 202 t 3 i 3 + t W2 t 323 + t i 2 t 323 — T 2 oiT 3 23^ 

-2 9 ^g 22 g 33 T 202 T 323 - 2s l yy 3 T 212 T323 + 2g 12 g' 2 g 33 T 212 T 323 . (A3) 
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